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Abstract 

The results of an experimental 
investigation on the combined effects of 
cooling channel aspect ratio and curvature for 
rocket engines are presented. Symmetrically 
heated tubes with average heat fluxes up to 1.7 
MW/m 2 were used. The coolant was gaseous 
nitrogen at an inlet temperature of 280 K (500 
°R) and inlet pressures up to l.OxlO 7 N/m 2 (1500 
psia). Two different tube geometries were 
tested: a straight, circular cross-section tube, 
and an aspect-ratio 10 cross-section tube with a 
45° bend. The circular tube results are 
compared to classical models from the 
literature as validation of the system. The 
curvature effect data from the curved aspect- 
ratio 10 tube compare favorably to the empirical 
equations available in the literature for low 
aspect ratio tubes. This latter result suggests 
that thermal stratification of the coolant due to 
diminished curvature effect mixing may not 
be an issue for high aspect-ratio cooling 
channels. 

Introduction 

The cooling passages of regeneratively 
cooled rocket engines have traditionally been 
designed with low height-to-w i dth 
dimensional aspect ratio. Recent experiments 
have suggested that improved thrust chamber 
liner cooling could be achieved by using high 
aspect ratio (AR>4) passages. 1 But, because 
only limited data were obtained in these 
experiments, further parametric studies were 
required to fully understand the potential 
benefits. 

The primary issue that remained from 
the previous work revolved around the 


curvature heat transfer enhancement. This 
enhancement provides additional cooling in 
the thrust chamber throat region, where the 
chamber liner receives the highest hot-gas side 
heat loads. The enhancement occurs because, 
as the coolant curves through the throat, a pair 
of secondary flow vortices develop which not 
only mix the coolant but also force cooler fluid 
to the hot wall. Since the previous experiments 
were performed with straight passages, the 
curvature effect was not investigated. 
Concerns have been raised that the tall thin 
passages might break down the vortex pair and 
reduce the curvature enhancement. A 
numerical investigation predicted that, if the 
secondary flow mixing were reduced, the 
coolant would become thermally stratified, and 
the hot-gas wall temperatures would be 
significantly higher. 2 

The issue was resolved in part by a flow 
visualization study of the secondary flow 
development in a curving high aspect ratio (AR 
= 5) channel. 3 The study showed qualitatively 
that the secondary vortex pair developed. But 
because the study was limited to low turbulent 
Reynold’s numbers and the strength of the 
secondary flow was not quantified, the results 
were not conclusive. 

In the current effort, a resistively 
heated tube experiment was initiated to obtain 
the detailed data necessary to model the 
curvature enhancement in high aspect ratio 
cooling passages. The coolant was gaseous 
nitrogen. Since the existing gas correlations 
are well founded for straight smooth-walled 
tubes, the investigation could be constrained to 
studying only the test section geometry effects. 
This report presents the results from tests 
which include validation of the experimental 
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facility with a circular test section and limited 
data from tests using a curved high aspect ratio 
(A R = 10) test section. The results from the 
circular test section are compared to the 
existing models for straight, smooth-walled 
tubes, and the curvature effects observed in the 
high aspect ratio test section are compared to the 
curvature factors offered in the literature. 


radius of curvature. Ito 8 later developed a 
similar curvature parameter for calculating 
the friction factor for turbulent flow in curved 
pipes 
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Theory 


The heat transfer coefficient in smooth 
straight tubes can, in general, be adequately 
predicted by models in the form of the Colburn 
equation. 4 Table 1 lists some of the variations 
that have been developed for straight smooth 
tubes with coolant properties evaluated at either 
the bulk or film fluid temperature. 4 » 5 
Simoneau and Hendricks 6 demonstrated that, 
for a gaseous coolant, the fluid properties could 
be grouped as a constant, and the heat transfer 
coefficient modeled by 

h = K{pV)l t d^ 2 ^ tJt: 

where K is a constant for a given gas, p is 
density, V is velocity, dh is the tube hydraulic 
diameter, and the subscripts b and w are bulk 
and wall conditions respectively. To apply 
any of these equations to a less than ideal 
situation, such as rocket engine design or 
analysis, several correction factors are 
usually applied such that 

Nu calc = Nu s V'Wr'Vc 

where the y’s are correction factors and the 
subscripts e, r, and c represent entrance, 
roughness, and curvature effects. The focus of 
this paper is \|/ c , the curvature factor, and the 
influence of cooling passage aspect ratio on it. 


where f c and f s are the curved tube and straight 
tube friction factors. 

Taylor 5 was one of the first to propose 
applying Ito’s 8 curvature factor to the heat 
transfer coefficient as 



The positive exponent is applied to the concave 
side of the tube and the negative exponent is 
applied to the convex side of the tube. 

The experiments of Niino, et al. 9 
revealed that although Taylor's factor is 
applied as a step at the entrance and exit of the 
curve, the actual curvature enhancement 
develops through the length of curve and 
persists beyond the exit of the curve. They 
obtained a better fit of their data with the 
following modification to Taylor's factor 


¥c = 


Re; 


\ 2 


i±o.02 r 




1 + tsin 


7T 


V 


L c + 15d h 


where x c is the distance from the start of the 
curve and L c is the length of the curve. Values 
for the coefficient b of 1 and 1/3 have been used 
to fit the data from two previous studies. 9 » 10 


A curvature parameter was developed 
by Dean 7 for laminar flow which was later 
successfully applied to predicting frictional 
pressure loss in a curving tube for laminar 
flow. Known as the Dean number, De, it is 
usually written as 


More recently, Kumakawa et al. 11 
proposed a slightly different version 



De = Re 



where Re is the Reynold's number and Rc is the 


The modeling capability of these 
factors will be compared to the experimental 
results. 
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Apparatus and Operating Procedure 

The flow system is shown 
schematically in figure 1. A vacuum jacketed 
run tank is located partially below the floor, 
and the coolant flows vertically through the 
system. The tank’s vacuum jacket allows for 
liquid nitrogen to be used as a coolant, but was 
not necessary for these tests because only 
gaseous nitrogen was used. The flow rate is 
redundantly measured by two venturies placed 
in series in the run tank dip tube, and bulk 
coolant conditions are obtained in mixing 
chambers immediately before and after the test 
section. The test section and mixing chambers 
are contained in a vacuum chamber to 
minimize external heat transfer. Heating is 
accomplished via the resistance of the tube with 
the electrical power provided by a direct current 
supply. 

Prior to a test, the run tank is charged to 
the desired pressure. When the test is initiated, 
the flow rate and back-pressure are controlled 
to set points by a programmable logic controller 
operating the respective valves. The power 
supply is then increased to the desired heat rate 
and the test section temperatures are allowed to 
stabilize. Once the temperatures are steady, the 
data are recorded, usually 10 cycles at 1 second 
intervals. The power is then cut off, and the 
tube is allowed to cool before stopping the flow to 
prepare for the next test condition. Further 
details of the test rig and operating procedures 
are described in reference 12. 

Test Hardware 

Figure 2 is a sketch of test section TSl 
and is representative of the instrumentation of 
test section TS5. Table 2 provides the 
dimensional details of both test sections. Test 
section TSl had a circular cross-section and 
was considered the checkout tube for validating 
the rig and data analysis procedures. Copper 
electrical busses were brazed to the tube for 
power connection, and extensive thermocouple 
and pressure tap instrumentation were 
applied. The thermocouples were spotwelded 
directly onto the tube, so it was necessary to 
take care in aligning the junction such that the 
test section voltage that is picked up by the 
thermocouple would be minimized. A 
description of the procedure used to correct this 


test section voltage induced error is presented 
in the analysis section. 

Test section TS5 had a cross-sectional 
aspect-ratio of 10 and a 45° curve in the plane of 
the larger dimension as depicted in figure 3. 
The curve radius to hydraulic radius ratio, R c 
/Rh , was 120, and the curve began 50 hydraulic 
diameters downstream of the first copper bus. 
This insured that the flow was fully developed 
from both a thermal and velocity standpoint at 
the entrance of the curve. Test section TS5 was 
fabricated by machining a curved open 
channel and electron beam welding the fourth 
side in place. One advantage of this technique 
over bending a straight tube is that the wall 
thickness is known, whereas it is not possible 
to predict exactly where the wall material will 
flow during the bending process. The wall 
thickness is critical in this type of experiment 
because it determines the local resistance and, 
hence, heat input to that portion of the tube, and 
ultimately the surface temperature. 

Analysis Procedure 

An analysis was conducted to calculate 
experimental convective heat transfer 
coefficients for the test sections from the 
thermocouple data. The analysis began by 
correcting the experimentally measured tube 
surface temperatures for errors due to the test 
section voltage picked up by the thermocouples. 
An initial series of tests were conducted at 
identical flow conditions but reversing the 
polarity of the voltage applied to the test section. 
Direct comparison of the temperatures 
indicated the extent of the error due to 
misalignment of the junction and the resulting 
potential picked up. Most of the thermocouples 
had negligible errors, but several required 
correction of the measured temperatures. No 
effort was made to measure the radiative or 
conductive losses from the thermocouple 
junctions, but a cursory analysis predicted 
only negligible losses. 

The corrected surface temperatures 
were then used to calculate inside wall 
temperatures assuming one dimensional 
conduction through the thickness of the wall. 
An iterative procedure that incorporates the 
variation of the thermal conductivity and 
electrical resistivity with temperature, similar 
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to that employed in reference 13, was used to 
obtain both inside wall temperatures and the 
local heat generation rate. Although the 
variation in heat generation rate along the 
length of the tube was determined from the 
measurements, the heat generation rate was 
assumed to be uniform through the wall 
thickness. The bulk temperature of the coolant 
at each station was calculated from the local 
heat input and the fluid properties obtained 
from the subroutine GASP. 14 A heat balance 
between the electrical power input and the 
coolant enthalpy gain was calculated as a 
check, to ensure that the heat losses from the 
system were not significant. For most of the 
tests, there was less than ± 3% discrepancy, but 
a few tests had a difference as large as ± 10%. 
It was believed the larger difference was due to 
the exit mixer not reaching thermal 
equilibrium, because of its large thermal 
mass, and not the result of a heat loss. Once it 
was determined that the heat balance was 
acceptable, the local experimental heat transfer 
coefficient and bulk property experimental 
Nusselt numbers were calculated. 

For test section TS5, the wall 
temperatures around the perimeter of the test 
section are not uniform, and the wall 
thicknesses (0.127 cm) are equal to the tube 
width. Therefore, a comprehensive analysis 
would incorporate three-dimensional 
conduction in the walls. However, in the 
present effort, the assumption of one- 
dimensional conduction was maintained. 
The impact that this assumption will have can 
be ascertained by considering the temperature 
gradients between the walls. In the straight 
thermal entrance region, the shorter walls 
were considerably hotter (on the order of 100 K) 
than the longer walls. Thus, heat will be 
conducted away from the shorter walls. 
However, the one-dimensional analysis 
assumes this conducted heat is convected to the 
coolant, and will calculate convective heat 
transfer coefficients that are slightly high. 
But, because the focus of this effort is the 
curvature effect, comparison of the curved 
section to the straight section should still be 
valid. 

Applying the assumption of one- 
dimensional conduction to the rectangular tube 
also required special treatment for the heat 


generated in the comers. This was handled by 
assuming that the heat generated in each 
corner was evenly distributed between the 
adjacent walls, which resulted in higher 
effective heat generation rates for the shorter 
walls, because the corner heat was distributed 
over a smaller area. This is in agreement 
with the hotter surface temperatures measured 
on the shorter walls. With these effective heat 
generation rates, the experimental heat 
transfer coefficients and bulk property Nusselt 
numbers could be calculated for each surface of 
test section TS5. 


Results and discussion 

The objectives of these test were to 
verify the current test and analytical 
procedures and to compare the predicted heat 
transfer coefficient curvature enhancement 
models available in the literature to the 
experimental enhancement in a high aspect 
ratio cooling channel. Test section TSl, a 
circular cross section tube, was used for the 
former and test section TS5, a curved, 
rectangular cross-section tube of aspect-ratio 
10, was used for the latter, A D-Optimal 
experimental test matrix was developed for the 
controlled experimental conditions: back- 
pressure, flow rate, and applied voltage. This 
D-Optimal design provides the benefits and 
efficiency of a designed experiment, but can be 
applied to an experiment where the design 
space is constrained. 15 In this case, the highest 
power input levels could not be operated with the 
lowest coolant flow rates without overheating 
the tube walls. The matrix also included 
several repeated tests. The range of test 
operating conditions for both test sections is 
shown in table 3. 

A representative wall temperature 
profile for test section TSl is plotted in figure 4. 
While the gradually increasing slope of the 
wall temperature in the central portion of the 
tube follows the increase in bulk temperature of 
the coolant, the drastic effects of the large 
copper power busses and the flow's thermal 
development are clear near the ends of the 43 
cm heated section. The heat transfer 
coefficients are also plotted in figure 4, and 
although, as expected, the end effects are 
strong, the behavior in the central portion of the 


4 



tube is typical of a gas-cooled tube. 

The raw data from twenty tests were 
reduced to obtain experimental Nusselt 
numbers using bulk fluid properties. The 
Nusselt numbers were then compared to those 
predicted by Taylor’s model in figure 5. 
Although the model included an entrance 
factor, only the central portion of the tube (from 
5 cm to 35 cm) was used for this comparison, to 
avoid obscuring the results with end effects due 
to the copper busses. The predicted Nusselt 
numbers are consistently higher than the 
experimental values, but fall within 20%. 

To improve the fit of the model to the 
data, a linear regression was performed on the 
wall-to-bulk temperature ratio exponent, and a 
value of -0.3 was obtained. The resulting 
model is compared to the data in figure 6 and 
shows excellent agreement. Examination of 
the data presented in reference 6 indicates that 
an exponent of -0.3 would fit the data better for 
wall-to-bulk temperature ratios less than two, 
as in this experiment. The ratio of 
experimental-to-calculated Nusselt number 
for the same model is plotted in figure 7 for all 
of the temperature stations, and shows the good 
fit in the central portion of the tube but not near 
the ends. 

Typical temperature profiles for the 
concave and convex walls of test section TS5 
are plotted in figure 8. For a straight tube, one 
would expect the wall temperature profile to 
follow the slope of the coolant bulk temperature 
profile, which is also plotted in figure 8. 
However, at the start of the curvature (~ X* = 11 
cm), the concave wall temperature turns 
dramatically downward. The convex wall 
profile becomes steeper at this point, but the 
change is not as dramatic. Unfortunately, the 
exit of the curve falls within the end effects of 
the copper power bus and so the wall 
temperature recovery could not be observed. 
However, the trends in both profiles are 
comparable to the curvature effects that have 
been observed in lower aspect ratio tubes. This 
indicates that strong secondary flows are 
present, and the coolant thermal stratification 
is not significant in the curve for these 
conditions. 

Finding an appropriate model for the 


curvature effects is more difficult than 
recognizing their presence. In figure 9 the 
ratio of experimental Nusselt number to 
calculated Nusselt number at the concave wall, 
including Taylor’s curvature factor, is plotted 
for five tests. Because the curvature 
enhancement increased through the curve, 
Taylor’s factor overcorrected at the entrance to 
the curve and undercorrected near the exit of 
the curve. However, the data all falls within ± 
20%. 

Figure 10 is a comparison of the 
experimental curvature factor from a 
representative test to that predicted by the 
methods of Taylor 5 , Niino, et al. 9 and 
Kumakawa, et al.n For the Niino et al. 
equation, the coefficient was assigned the 
value 1/3. Both the concave and convex surface 
curvature factors are shown. Note that the 
predicted values for the two surfaces are not 
symmetrical because the concave side of the 
curve is longer than the convex side. Although 
both the method Taylor and Kumakawa, et al. 
underpredict the experimentally derived 
values for the concave surface, Kumakawa’s 
method models the development much better. 
The method of Niino, et al. over predicted the 
curvature factor, even though the smaller value 
of the coefficient b was used. No attempt was 
made to obtain a better fit of the data, because as 
noted earlier, the one-dimensional analysis 
used to determine the experimental Nusselt 
numbers may have introduced some errors. 
Therefore, it was not determined whether the 
discrepancies were real or a result of the 
simplified analysis. The curvature effect on 
the convex surface is overpredicted by all three 
models. 

Concluding Remarks 

An electrically heated tube facility was 
used to investigate the combined effects of 
cooling channel aspect ratio and curvature. 
For all of the data presented, gaseous nitrogen 
was the coolant. 

To validate the facility, a single 
circular cross-section tube was tested, and the 
results compared favorably with the classical 
heat transfer correlations available in the 
literature. 
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The results from a curved rectangular 
cross-section tube of aspect ratio 10 showed that 
a strong curvature induced heat transfer 
enhancement was present even at this high 
aspect ratio. The implication of this result was 
that secondary flows developed that minimized 
the coolant thermal stratification. The data 
also indicated that Taylor’s curvature factor, 
\|r c , underpredicted the experimental results by 
up to 20%. However, the one-dimensional 
analysis that was applied to the data may have 
calculated a higher experimental heat transfer 
coefficient than actually occurred and caused 
the difference. 

The development of an optimum factor 
for accurately predicting the curvature heat 
transfer enhancement in high aspect ratio 
channels will require further testing and more 
comprehensive analysis of the data, as the 
results presented here were limited to a single 
curvature and a single aspect ratio. However, 
based on these results, the author recommends 
that, for coolants for which the Colburn type 
Nusselt number models apply, and up to an 
aspect ratio of 10, Taylor’s factor, or the 
modified version suggested by Kumakawa, 
should be used to predict the curvature 
enhancement in cooling channels. 
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Table 1. Heat Transfer Coefficient Models For Straight, Smooth-Walled Tubes 


Dittus - Boelter 4 

Nu s = 0.023 Re* 8 Pr b 04 

Taylor^ 

Nu s = 0.023 Re° 8 Pr° 8 

( T ^ 

* w 

InJ 

-f 0.57- 159 1 
l X/D) 

Sieder-Tate 4 

Nu s = 0.027 Re° 8 Pr° 8 

' 

to 

0.14 


Table 2. Test Section Details 


Test Section 
Number 

Material 

Inside 

Dimension (cm) 

Outside 

Dimension (cm) 

Heated 
Length (cm) 

Curvature 

TSl 

SS 304 

0.564 diameter 

0.635 diameter 

43.2 

None 

TS5 

Inconel 718 

0.127 x 1.270 

0.381 x 1.537 

25.4 

45° Bend with 
13.7 cm mean 
radius 


Table 3. Ranges Of Operating Conditions For Each Test Section 


Parameter 

Test Section TSl 

Test Section TS5 

Bulk Reynolds 
Number, Re b 

400,000 to 1,500,000 

250,000 to 900,000 

Wall-to-Bulk 
Temperature Ratio, 
T w /T b 

lto2 

lto2 

Back Pressure 
(MPa) 

2.1, 4.8, 8.3 

2.1, 4.5, 6.9 

Coolant Bulk Inlet 
Temperature (K) 

-280 

-280 

Average Heat Flux, 
q” ( MW/m 2 ) 

0.30, 0.84, 1.70 

0.34, 0.74, 1.31 
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Figure 1. Heated tube facility schematic. 


Figure 2. Test section TS 1 instrumentation. 



Figure 3. Test section TS5 geometry. 



Figure 4. Typical experimental wall temperature (T ), bulk temperature (T^'j, 
and heat transfer coefficient (h ) profiles for test section TSl. Bulk temperature 
is adjusted by a constant 100 K for the plot. 
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Figure 5. Comparison of Nusselt numbers calculated by Taylors equation for straight 
smooth-walled tubes to the experimental values for test section TS1 cooled with gaseous 
nitrogen. 



Figure 6. Comparison of Nusselt number calculated with a modified version of 
Taylor's equation to the experimental values for test section TSl. 
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Figure 7. Nusselt number ratios plotted for all temperature stations of test section 
TSl. 



Figure 8. Comparison of typical wall temperature profiles for the concave 
and convex sides of test section TS5 to the coolant bulk temperature. 
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Figure 9. Nusselt ratios for the concave wall of test section TS5 using 
Taylor’s curvature factor. 



Xj, inches 

Figure 10. Comparison of curvature factor calculated by the models to typical 
experimental results. 
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